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SUMMARY 

A numerical study for confined, axlsymmetrlcal , turbulent diffusion flames 
Is presented. Local mean gas properties are predicted by solving the appropri- 
ate conservation equations In the finite-difference form with the corresponding 
boundary conditions. The k-e two-equation turbulence model Is employed to 
describe the turbulent nature of the flow. A two-step kinetic model Is assumed 
to govern the reaction mechanism. The finite reaction rate Is the smaller of 
an Arrhenius type of reaction rate and a modified version of eddy-breakup 
model. Reasonable agreement Is observed between calculations and measurements, 
but to obtain better agreement, more work Is needed on Improvements of the 
above mathematical models. However, the present numerical study offers an 
Improvement In the analysis and design of the gas turbine combustors. 


INTRODUCTION 

The design of modern gas turbine combustors Is an extremely complex proc- 
ess. Many factors can Influence their performance, Including combustor geome- 
try, flow properties of oxidant and fuel, flow turbulence, and combustion 
process. Consequently, methods that could aid In understanding and predicting 
these phenomena would be very useful In designing more efficient combustors. 
This report describes numerical calculations of a turbulent reacting flow In a 
combustion chamber In which the flow Is confined and axlsymmetrlc with sudden 
expansion. 

The reacting flow fields considered are gaseous-fueled and are governed 
by a set of time-averaged conservation equations. Air and natural gas are used 
as the oxidant and fuel, respectively. In order to close these governing equa- 
tions, a set of mathematical models are required for the description of various 
physical processes. Neglecting the radiation effect, two main features need to 
be modeled: flow turbulence characteristics and chemical reaction. The k-e 

two-equation turbulence model Is employed for the description of eddy proper- 
ties. For the chemical reaction rate, a modified version of eddy-breakup model 
Including kinetic effect, Is used. 

The objective of present study Is to predict the chemical compositions and 
to compare with measurements and then examine the validity of the mathematical 
models for combustor problems. It Is Impossible to obtain analytical solu- 
tions, so the computational fluid dynamic techniques are required. A TEACH 
computer code using the Bounded Skew Hybrid Differencing (BSHD) Is used In the 
present study (refs. 1 and 2). The Pressure-Implicit Split Operation (PISO) 
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predictor-corrector procedure (ref. 3) is employed for the calculation of pres- 
sure distribution. 


GOVERNING EQUATIONS AND BOUNDARY CONDITIONS 

The partial differential equations governing the flow considered herein 
are of the form: 

[fe (oW * t!f ’■<'>''♦>] ■ ( b It)* r If ( rb I?) ♦ s * (1 > 

where the corresponding values of dependent variable <t>, effective exchange 
coefficient b, and source term S^ are given In table I. The equations 
represent conservation of mass, three components of momentum, turbulent kinetic 
energy, turbulent dissipation rate, total enthalpy, mass fraction of fuel, mass 
fraction of CO, and the mixture fraction. The mixture fraction f which Is 
defined as 


can be considered as the total mass fraction of burned and unburned fuel. r 3 
Is the mass stoichiometric ratio of H 2 O and fuel In the chemical reaction 
equation. 

The effective viscosity Is given by 

v e ff - n + vi 

where an< * w are the eddy and laminar viscosities, respectively, vt 
Is calculated from the solution of the k-e turbulence model which will be 
discussed later. The turbulent Prandtl/Schmldt numbers, °fu» °co» an( * °f» 
are assumed to have a value of 0.9. The source terms, Sf u , S co , and r 2 , In the 
mf u and m c0 equations are Involved with the combustion model and will be 
discussed In the later sections. 

Since the equations are elliptic In form, boundary conditions are required 
at all the boundaries. Figure 1 shows the geometry of the combustor consid- 
ered. The flow measurements for this geometry were made by Lewis and Smoot 
(ref. 4). Some boundary conditions can be obtained from their measurements. 
Table II shows the measured boundary conditions and the compositions of air and 

fuel. The Inlet velocities are parallel and assumed to have power law distri- 
butions. This Implies that the flow Is nonswirling and the tangential momentum 

equation can be neglected. The Inlet turbulent kinetic energy for air flow was 

given as 11.765 m^/s? and the corresponding value for fuel flow was meas- 
ured to be 1.633 m 2 /s*. The dissipation rates are estimated by the follow- 
ing relation: 


e = 
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where Is the diameter of the flow passage. At all walls, the nonslip 
condition Is applied and temperature Is provided by the measurements. The 
gradients of other dependent variables are assumed zero. The dependent vari- 
ables at the wall are linked to those at the grid node next to the wall by the 
logarithmic law of the wall. The symmetry axis provides a boundary condition 
of the form 3<|>/ar = 0 except the radial velocity v which Is zero at the 
centerline. At the exit, all gradients are assumed to be zero, l.e., 

3<|>/dx = 0. 


THERMODYNAMIC PROPERTIES 

The local density of gas mixture Is obtained from the equation of perfect 
gas 


P = 



where M Is the mixture molecular weight calculated from the relation 



where m^ and are the mass fraction and molecular weight, respec- 
tively, for the 1th species. The total enthalpy H Is defined as 
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usually neglected 
for low-speed flow 


where Hf u and H c0 are the heats of reactions for fuel and CO, respectively, 
corresponding to the two-step overall kinetic equations, and the species spe- 
cific heat Is calculated by 


= a 1 + b^T + c^ 2 + d J| T 3 


(3) 


The constants a^, b^, c^, and d^ can be found from reference 5. The temper- 
ature T Is calculated Iteratively by using Eqs. (2) and (3). 


CHEMICAL KINETICS MECHANISM 

The simplest chemical mechanism Is the one-step overall kinetic mechanism 
which Is the most convenient for numerical modeling. This approach assumes the 
oxidation process to occur directly to CO 2 and H 2 O. However, It does not 
account for the characteristics of hydrocarbon oxidation since small amounts 
of Intermediate hydrocarbons and somewhat larger amounts of CO are usually to 
form prior to significant productions of CO 2 and H 2 O. 
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The next stage of complexity Introduced has been the two-step overall 
kinetic mechanism which separates the highly exothermic oxidation of CO to 
CO 2 from the less exothermic oxidation of the hydrocarbon to CO (refs. 6 
and 7): 


C x H y + 


(f + ?)(° 2 + nN 2> ^ XC0 + 2 H 2° + (2 + 4) 


nN, 


XCO ♦ | (0 2 + nN 2 ) -» XC0 2 + | nN 2 . 


In the first reaction, 

ri = (mass of 0 2 )/(mass of fuel), 
r 2 = (mass of C0)/(mass of fuel), 
r 3 = (mass of H 2 0)/(mass of fuel); 

In the second reaction, 

r 4 = (mass of 0 2 )/(mass of CO), 
rs = (mass of C0 2 )/(mass of CO). 

The values of these ratios can be obtained In the following: 



Here the W's are the molecular weights of the chemical species. With the 
values of r's, we can calculate the mass fractions of 0 2 , C0 2 , and H 2 0 from 
the following: 

"*0 2 = R(1 ' f) + r 4 m C0 + (r l + r 2 r 4 )(m fu " R = °* 233 
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It Is noted that the fuel considered In the present study Is a natural gas con- 
sisting of several components such as CH 4 , C 2 H 6 , etc., so the calculations of 
r's described above need to be modified. This can be done without much 
difficulty. 


TURBULENCE AND COMBUSTION MODELS 

The k-e two-equation turbulence model Is used for the closure of 
Reynolds equations. The eddy viscosity Is obtained from the relation 

C P k 2 

.. _ -y — 


Values of the constants associated with the calculation of eddy viscosity using 
the k-c turbulence model are given In the following (ref. 8): 

C 1 = 1.44, C 2 = 1.92, C^ = 0.09, K = 0.42, o k =1.0, 



In the near-wall region, equations were Introduced to link the values of 
dependent variables on the wall to those In the logarithmic region (ref. 9). 

This turbulence model has been tested previously for many recirculating 
flows (refs. 10 to 13). The Investigation Indicated that acceptable agreement 
was obtained between measured and predicted flow pattern, flame structure, and 
heat flux. 

Equations are required for the reaction rate of fuel and CO. Several 
models have been proposed to predict burning rate of fuel In turbulent environ- 
ments. In the present study, an Arrhenius type reaction rate Is compared with 
a modified eddy-breakup rate, and the smaller of the two rates controls the 
reaction (ref. 7). The corresponding rates for the first reaction Is expressed 
as: 


where 


Sfu = -Min. (S-j ,S 2 ) 
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S-| Is the Arrhenius rate of fuel oxidation as controlled by chemical kinetics. 
The eddy-breakup model of Spalding expresses the rate of oxidation as Influ- 
enced by turbulence Intensity and scale, and concentration of unburned fuel. 
This model Is applicable to premixed flames. However, the combustion In gas- 
turbine combustors Is neither fully premixed nor entirely diffusion controlled, 
so a term Is added In the model to determine the rate of fuel oxidation as con- 
trolled by the availability of the oxygen. 

Similarly, the reaction rate of CO can be expressed as: 

S^q = -Min. (S3, S4) 

where 



1-) and I2 are the stoichiometric ratios of the first and second 
chemical equations, respectively. The eddy-breakup constants, Cr 1 and 
CR t 2 are set to be 2 *° and respectively. The other constant appearing 
In’the Arrhenius equations are the following values (refs. 14 and 15): 


F, = 3.3X10 14 , — = 27 000 K, 

1 R 


F, = 2.2xl0 8 , — = 12 500 K, 

£ R 


SOLUTION PROCEDURE 

The governing equations and the associated boundary conditions are solved 
by a TEACH computer code using the Bounded Skew Hybrid Differencing (BSHD) 
method (refs. 1 and 2). The pressure distribution Is estimated by the 
Pressure-Implicit Split Operation (PIS0) predictor-corrector technique 
(ref. 3). Calculations are performed with a nonuniform grid distribution com- 
prising 60 by 47 nodes with concentration of the nodes In the near-wall region 
and centerline regions. The values of the underrelaxatlon parameters for each 
dependent variable $ are listed In table III. The solutions are assumed 
to be convergent when the maximum residual Is reduced by three orders of magni- 
tude for any node for any <t> equation. The numerical calculations are per- 
formed on a CRAY1 computer with CPU time of around 56.7 min. About 3960 Itera- 
tions are required to obtain convergent solutions for the present calculations. 
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RESULTS AND DISCUSSION 


The numerical predictions compared with measured data are shown In 
figures 2 to 9. From figure 2, we observe that the calculations and measure- 
ments for the mixture fraction distributions are In good agreement except In 
the regions near the centerline and In the recirculating regions close to the 
corner. The calculations overpredict and underpredict, respectively, the mix- 
ture fraction In these two regions. However, It Is noted that In the develop- 
ing region near the centerline, large radial concentration gradients exist. 

The measured (or collected) sample represents the average species concentration 
over the area of the probe Inlet, so the magnitude of the error can be up to 
22 percent while the error Is negligible for the radial positions outside 3 cm 
or for axial positions aft 25 cm (ref. 4). With these measurements corrected, 
the discrepancies between calculations and measurements are expected to be 
reduced In the developing region. The calculated CH 4 and O 2 concentration 
profiles In figures 3 and 4, respectively, show that the reaction rates are 
slower than observed, especially In the two regions mentioned above. Conse- 
quently, the concentration of CO 2 and H 2 O (figs. 5 to 8) are underpredicted 
In both the recirculating region near the corner and developing region In the 
centerline. In far downstream, the predicted H 2 O profiles are somewhat 
higher than the measured data probably due to the neglect of H 2 as one of the 
Intermediate products. The concentration profiles of CO 2 and CO (figs. 7 
to 9) Indicate that the reaction rate for CO 2 Is too low and that the radial 
spread of CO profiles Is underpredicted. Two temperature profiles are shown 
In figure 9. Since the reaction model used In the present study underpredicts 
the chemical reaction rate, the calculated temperature Is lower than the meas- 
ured data In the centerline developing region. 

According to the above Investigations, we can conclude that the overall 
discrepancy between the calculations and the measurement are due to the Inap- 
propriate mathematical modeling In the centerline developing region and the 
corner recirculation region. Except for these regions, acceptable agreement 
between numerical and measured data Is made using the k-c turbulence model 
and reaction model used In this study. 


CONCLUSIONS 

In the present study, a numerical calculation for turbulent combusting 
flow In a confined sudden expansion gas-turbine combustor has been Investi- 
gated. The local mean gas properties are predicted and compared with measure- 
ments. The species concentration distributions are reasonably well predicted 
except In the centerline developing region and corner recirculation region. 
Some factors which may cause prediction error are: 

1. Inappropriate combustion model. 

2. Weakness of turbulence model. 

3. Inappropriate chemical mechanism by using two-step overall kinetics. 

4. Neglect of minor Intermediate species. 
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5. Assumption of equal dlffuslvltles of all chemical species. 

6. The numerical error from the computer code Itself. 

It Is believed that the major discrepancy between calculations and measurement 
Is due to Items 1 and 2. Research In these areas Is still a challenge to us. 
However, more accurate turbulence model and combustion model are an urgent need 
In Improvement of numerical study. 
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TABLE I. - GOVERNING EQUATIONS CORRESPONDING TO EQUATION (1) 



TABLE II. - COMBUST OR PARAMETERS AND BOUNDARY CONDITIONS 
Primary gas: naturaT gas 


Temperature, K 300 

Flow rate, g/s 2.84 

VeTocity , m/s 21.3 

Composition, mol % 

CH 4 83.9 

C 2 H 6 7.1 

N 2 2.4 

C0 2 1.3 

Ar 5.3 

Secondary gas: air 

Temperature, K 589 

Flow rate, g/s . 36.3 

Velocity, m/s 34.3 

Composition, mol % 

N 2 79 

0 2 21 

Pressure, N/m 2 94 000 

Wall temperature, K 1140 




TABLE III. - VALUES OF UNDERRELAXATION PARAMETERS 


Variable 

Underrelaxatlon 

parameter 

U velocity 

0.3 

V velocity 

.3 

Primary pressure correction 

1 .0 

Turbulence energy 

.4 

Dissipation rate 

.4 

Total enthalpy 

.6 

Fuel concentration 

.6 

CO concentration 

.6 

Mixture fraction 

.6 

Second pressure correction 

.7 

Viscosity 

1 .0 

Density 

.1 



r 1 - 0.008 M 

r 2 = 0.0111 n 

r 3 = 0.0286 M 
R = 0.1016 M 
L = 1 ,52*1 n 

FIGURE 1. - GEOMETRY OF AXISYMMETRIC COMBUSTOR WITH COAXIAL FUEL AND 
AIR JETS. 
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